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THE STRAIN OF A GRAVITATING, COMPRESSIBLE ELASTIC 

SPHERE* 

BY 

L. M. HOSKINS 

In the solution of the problem of the strain of a homogeneous elastic sphere, 
as given by LamS and by Kelvin, the bodily forces are supposed expressible 
as known functions of the coordinates of position. When self-gravitation is 
among the forces considered, however, the bodily forces depend in part upon the 
state of strain, and it is only in the case of assumed incompressibility that 
the usual method of solution is applicable. 

The present paper has for its main object the solution of the problem of the 
strain of an isotropic elastic sphere, initially homogeneous, due to disturbing 
forces of a certain type, taking into account the changes in the gravitational 
forces which result from the strain. 

The contents of the paper fall under the following four heads : 

I. Strain of a gravitating, compressible elastic sphere under the action of 
tidal or centrifugal forces, with numerical computations for the case in which 
the surface is free from stress. 

II. Strain of a gravitating, compressible elastic sphere covered by a shallow 
ocean, under the action of tidal or centrifugal forces. 

III. Effect of compressibility upon estimates of the rigidity of the earth. 

IV. Strain of a gravitating, compressible elastic sphere under the action 
of small disturbing forces derivable from a potential which is any spherical 
harmonic of degree not less than 2. 

I. Strain of a compressible, gravitating elastic sphere under the 

ACTION OF tidal OR CENTRIFUGAL FORCES. 

§ 1. Physical analysis of the problem. 

The object of the analysis is to determine the strain produced in an elastic 

sphere by the action of disturbing forces of the type of tidal or centrifugal forces 

— i. e., forces derivable from a potential proportional to P^i"^ , in which P^ is the 

zonal surface harmonic | cos^ 6 — | . The disturbing forces are assumed small 

in comparison with the actual forces of gravity, but the changes in the gravita- 

* Presented to the American Mathematical Society (San Francisco), September 25, 1909. 
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204 L. M. HOSKINS: THE STRAIN OF [April 

tional forces due to the strain are assumed to be of tlie same order as the dis- 
turbing forces. This accords with the facts in the case of the actual tidal and 
centrifugal forces acting upon the earth. 

Before the action of the disturbing forces the body is assumed to be in 
equilibrium under the action of the mutual gravitation of its parts without sur- 
face stress. The stress-condition at any interior point need not be specified 
more definitely ; it may be hydrostatic equilibrium, or any departure from that 
condition which is consistent with the forces of self-attraction. The stresses 
and strains considered in the following analysis are to be regarded as super- 
posed upon the supposed initial condition, it being assumed that small depar- 
tures from this initial condition follow the laws of elasticity as ordinarily given 
for a homogeneous isotropic solid.* 

When the external disturbing forces come into action, the body assumes a 
new configuration of equilibrium. The gravitational force which now acts upon 
a given element of volume is different from that which originally acted upon the 
same element, for two reasons, — because the density of the attracted element is 
changed, and because the density-distribution of the attracting mass is changed. 
In computing the new value of the gravitational force the new configuration of 
the body may be treated as made up of two parts, — the original configuration, 
and a certain mass-distribution superposed upon it. The total attraction upon 
a volume-element may then be expressed as the resultant of the following parts : 

(1) The attraction of the original mass-configuration upon the original mass of 
the element ; (2) the attraction of the original mass-configuration upon the incre- 
ment of mass of the element ; (3) the attraction of the increment of mass-distri- 
bution upon the original mass of the element ; and (4) the attraction of the incre- 
ment of mass-distribution upon the increment of mass of the element. Of these 
four forces (1) is balanced by the initial stresses and may be omitted from the 
equations of equilibrium of the element ; (2) and (3) are of the same order of 
magnitude as the disturbing forces ; while (4) is negligible in comparison with 

(2) and (3). Hence (2) and (3) are the only bodily forces to be considered in 
addition to the external disturbance. 



* It has been held by certain writers that the existence of initial stresses of great magnitude 
vitiates the application of the theory of elasticity to a gravitating body comparable in size and 
mass with the earth. As to applications which are sometimes made of the theory, the objection 
appears to be undoubtedly valid. For example, it would be absurd to compute the actual 
stresses within the earth as it they were due to gravitation acting upon a body initially free 
from stress ; even if such an assumption were in harmony with any reasonable supposition as 
to the history of the earth, the strains thus computed would far surpass the limits of elasticity. 
This objection does not, however, appear to hold against such an application of the laws of 
elasticity as is made in the present paper. To assume that a body under great stress, whether 
this is everywhere purely a normal pressure or whether involving also tangential stresses, is 
elastic for small departures from the initial configuration, is not inconsistent with what is known 
regarding bodies under ordinary experimental conditions. It can hardly be doubted that the 
great body of the earth is elastic as regards slight changes of configuration. 
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The original density being called p, the increment of density Sp is equal to 
— pA, if A is the cubical expansion. Hence the force (2) is to be expressed as 
the attraction of the undisturbed mass upon a density — pA, while (3) is the 
attraction of a mass-distribution of density — pA upon a density p , together with 
the attraction due to the surface inequality produced by the strain. 

Besides the bodily forces it is necessary to consider the stress-condition at the 
surface. Upon the actual bounding surface it is assumed that the stress remains 
zero ; but if (as is convenient) the equations of interior equilibrium be assumed 
to apply throughout the original spherical volume, the surface of this volume 
must be regarded as sustaining a stress whose value must be properly expressed. 
It is shown below that, to the first order of small quantities, this stress is radial 
and equal to the weight of the surface inequality produced by the strain. 

§ 2. Formation of equations of equilibrium for an element of volume. 

(a) General form of equations in terms of polar coordinates. Let the axis 
of symmetry of the disturbing forces be taken as initial line of a system of polar 
coordinates, the displacements in the three primary directions at a point 
{r, 6, (fi) being called 

U^ , Ug , u^ , 

while A is the cubical expansion. In the case now under consideration 

M, = 0. 

Let the total bodily force per unit volume have components 

(the </)-component being 0). 

In accordance with a common notation, let fi denote the modulus of distortion 
and A. + |/U. the bulk-modulus of the material. The equations of equilibrium 
for an elementary volume then assume the following forms : * 

(1) iX + 2,)-^^--r^^^^^--^)^e\+pF = ^, 

(2) (X+2M),;^+-{b^4(^,/-^---.-g^)sin^J + ,i., = 0; 

in which 

1 r 5 d "1 

(3) ^ = -2-gi„-^ [ Q^ ( »•' w. sin ) -f -^^ ( rw, sin ) J . 



(6) Assumptions. It is found that by making certain assumptions the equa- 
tions are reduced to a form admitting of solution. It will be well to state these 

*See A. E. H. Love's A Treatise on the Mathematical Theory of Elasticity, 2d edition (1906), 
pp. 138, 56. The terms representing bodily forces are omitted from the equations as given on 
p. 138. 
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assumptions explicitly, at the same time remarking that their justification is 
the fact that they lead to a solution satisfying all the conditions of the problem. 

(1) The potential of the disturbing forces being proportional to r'^P^i in 
which P^ = I- eos^ d — ^^ , it is assumed that the radial displacement and the 
cubical expansion depend upon 6 in the same way as this potential does ; i. e., 
each is assumed to be the product of P^ by a function of r. 

(2) It is assumed that the mean radius of an originally spherical surface is 
unchanged by the strain. 

(c) Linear displacements replaced hy ellipticity and maximum angular dis- 
placement. It is found convenient to replace u^ and u^ in the foregoing equa- 
tions by two quantities e, a, defined as follows : 

The former represents the ellipticity of the surface originally spherical and 
of radius r ; that is, 

(4) u^ = - lerP^ 

(it being assumed that a positive value of e corresponds to oblateness). 

The quantity a represents the angular displacement of a radius vector for 
which = 45°. It is then found that, for any point (r, ^), 

(5) '^ = a sin 2^. 

This follows from the assumption that A involves Q only in a factor P^. 

To make the last statement clear, notice that the introduction of e in place of 
u^ reduces (3) to the form 

(^ <le,\ ^ 1 5(M«sin0) 

\ ^ dr J ^ r sin t' oB 

In order that the last term may be equivalent to P^ multiplied into a function 
of ?•, it is easy to show that we must have 

• on 

u„ = a sin W + -; — Ti , 
* sm ^ ' 

in which a and /3 are independent of 0. Infinite values of Ug being excluded, 
/3 must be .* 

If, now, u^ and u^ be replaced by e and a by means of (4) and (5), equations 
(1) and (2) take the following forms : 

*It is otlierwise evident that, tor a given value ot r, «# will have equal values for comple- 
mentary values of 6 . This is seen by considering the effect of superposing strains of the same 
type about different axes. Thus the ^-displacements in a given meridian plane (say that of xz) 
will be annulled by superposing upon the given strain either (a) an exaotly similar strain sym- 
metrical about the r-axis, or (6) a strain symmetrical about the z-axis which is just the reverse 
of the given strain. 
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(7)[/*''$+4M^-2(3\+5M)^+(X+M)J+(3X+5M)^]sin20+pi^,=O, 

while the value of A becomes 

(8) A = (4«-2e-f4')A- 

(d) Expression of bodily forces. The bodily forces acting upon an element- 
ary volume are (1) the disturbing force (here supposed to be the centrifugal 
force due to a uniform rotation about a fixed axis), (2) the gravitational force 
due to the change of density of the attracted element, (3) the attraction due to 
the change of configuration of the attracting mass. We proceed to get expres- 
sions for these in terms of r , 6 , e and a . 

(1) The disturbing force. 

Representing the angular velocity of the rotation by o), the centrifugal forces 
have a potential 

(9) W iwV^P^ 

(neglecting a term independent of 6). If a is the radius of the sphere, g the 
value of gravity at the surface, and c the ratio of equatorial centrifugal force to 
gravity, i. e., 

(10) c_'^°, 
we may write 

(11) ^=-fa^P^- 

The parts of F^ and F^ corresponding to this potential are 

(2) Force due to change of density of attracted element. 

The attraction per unit mass at any point, due to the original homogeneous 
sphere, is a radial force equal to 

r 

The increment of density of the attracted element being — pA , the force per 
unit volume is 

/srA 
a 



(-/-)(-, A) = ,' 
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hence this force is accounted for by adding to the value of pF^ a term 

(3) Force due to change of configuration of the attracting mass. 

We have now to express the value of the potential at any point due to a 
mass-distribution of density — pA throughout the original spherical volume, 
together with that due to the surface layer superposed upon this sphere by rea- 
son of the strain. 

Consider a shell of radius r', thickness dr, and density 

The potential at a point (r, 6) due to this shell is 

/ de \ r^ 

(15) _4^T,p(^4a'-2e'-|r'^,j^cZr'-P, (if»-<r'), 

(16) _|^^p(^4a'-2e'-fr'j;)^c^/.p, (ifr>r'). 

Hence the potential at any point ( r , ) within the spherical volume of radius a , 
due to the density-distribution — /sA , is 

\-7T^pP^ r 1 



(17) 



{^'/(^'^'-^^'-t^'S^)^"^^' 



^.f(--^^'-t^'£)?]- 



To this must be added the potential due to the surface inequality. If c, is the 
surface value of e , the thickness of the layer is — |ej aP^ , and the potential due 
to it at a point (»•, ^) within the body is 

(18) _J_^^/,e,r^P,. 

The total potential due to the change of mass-distribution may therefore be 
written in the following form, in which g replaces ^Trypa, the surface value of 
gravity : 

3^ A I 



(19) 



5a 



{^jr'(4,'-2.'-,4,)/'*' 



f" / de \ dr' ~1 

By partial integration we find 

1 r' / de \ 4 /" 
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so that the expression for V reduces to 

The forces derived from the potential Fare (per unit mass) 

<^>) ^?-'^'['.r(''--+«>"*--x>«-«'>7-8-]- 

(4) Components of total bodily force. 

Combining the foregoing results, we have for the components of the bodily 
force per unit volume at the point (r, 0), 

(23) 

+ gJ^f3«'4-e>'^c?/-^-^J"(2a'-0^^], 

(24),i.,=^^^[.c.+ |,J>«'4-e>'^..'4^f(2a'-e')^]. 

(e) Final form of equations of equilibrium for elementary volume. In 
writing the two equations of equilibrium for a volume element in their final form, 
it is convenient to introduce the factor r throughout each equation, and then to 
replace r as independent variable by x = rfa. The equations then take the 
following forms : 

ft (Y fi P fiP 

(25) + pga — icx^ + f 4a — 4e — |x ,^ ja;^ 

da GiOi Q/B 

'^'^'d^^'^^'^dx -2(3X + 5M)a + (X + /.)x^~^ + (3X + 5M)e 
(26) 

+ pgaYlcx' + ^j {Za+e)x'dx+^ j (2«'- e') ^ J =0. 
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§ 3. Surface conditions. 

Upon the bounding surface of the strained body the stress is zero. It seems 
most convenient, however, to treat the strained body as extending to or beyond 

the original spherical boundary r = a'va. every direction, 
and to determine what stress must exist upon this sur- 
face in order that the surface r ■=■ a{\ — ^e^P^) may 
be free from stress. For this purpose consider the 
layer of matter of density p included between the sur- 
faces r = a, r = ai'l — §6,^2). Let equations of 
equilibrium be formed for an elementary portion of this 
layer bounded by a meridian plane and a plane parallel 
to it at an infinitesimal distance h , and by two radial 
planes perpendicular to the meridian inclined to each other at an infinitesimal 
angle i|r ; the section of such an element by the meridian plane is represented by 
ABB' A (Fig. 1). 
We have 

AB = af, AA'=-iae,B„ 




Fig. 1. 



d 
BB = AA' + ■f j^{— %ae^Pi) = AA + aCji/r sin 



2^. 



The stresses per unit area at the point A are 

(27) f>- = XA + 2M^'^=[4x«-2(\ + |^)e-f(X + 2/.)a5^~]p„ 

= 4\a-2(X+ %ii)e— i'^x-^ Pj + 4,itacos2^, 



(28) 



(29) 



->, /^Wfl W(i 1 d)i \ ( da\ . „^ 



From these expressions values may be found for the total normal and tangential 
forces on the faces AB, AA., BE of the element, while the total bodily force 
is the weight of the element, which is equal to pg into its volume. By writing 
the two equations of equilibrium for the element, and retaining only terms of 
lowest order in the small quantities e, a, it is found that 

(30) -fr+pg{\eaP,) = ^, 

(31) rO = 0. 

The equations to be satisfied at the surface r = a are therefore 
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(32) -p^ = 4Xa-2(X + f^)e-|(\+2M)x^- = fP5'ae, 

<^«) si^.-(' + »£) = "■ 

§ 4. Preliminary discussion of special cases. 

The problem is now reduced to that of solving the simultaneous ordinary 
differential equations (25) and (26), subject to the conditions expressed by (32) 
and (33). Before proceeding to the general solution, it is of interest to note 
the form assumed by the equations in special cases for which the solutions are 
known. Such cases are the following : 

(A) A non-gravitating body acted upon by bodily disturbing forces having a 
potential proportional to Pjr* (or free from bodily forces), and 

(1) With any given surface stresses of type corresponding to the bodily 
forces (i. e., the radial stress proportional to f^ and the tangential stress to 
sin 26), or 

(2) With any given surface displacements of the same type (i. e., the radial 
displacement proportional to P^ and the tangential to sin 26). 

[B) A gravitating incompressible body with bodily and surface forces and 
surface displacements as in any of the cases under {A). 

Of these cases the complete solution will here be given only for that of a non- 
gravitating body with surface free from stress. The modifications in the 
analysis required by the other cases under (A) are sufficiently obvious, while 
the corresponding case of (S) is given later on as a limiting case of the general 
solution for a gravitating compressible body. 

It may be remarked, also, that after the general solution has been given it 
will be obvious how it must be modified if the surface conditions are changed. 

(a) Solution for the case of a non-gravitating body acted upon by bodily 
forces but free from surface forces. In this case all the terms in equations 
(25) and (26) which represent bodily forces disappear except those representing 
the disturbing forces, and the equations become 

(34) 4(X+M)x£ -8/.a-|(X + 2M)x^^-f(X+2/.)x^ + i/^e-ygacx'^O, 

(35) /^a''^ +4/^»'^ -2(S\+5fj,)a+(K+fj,)x~+{S\+5iJi)e + ^pgacx" = 0. 



The surface conditions (32) and (33) become 

di 
' dx~ ^ '' "^ "'" "^ dx 



(36) 4Xa-2(\-j-|yn)e-|(X + 2;u.)x~ = 0, e + X;^ = (a; = l). 
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It is easily shown that these equations may all be satisfied by values of e and 
a of the forms 

(37) e = A^ + A.^x\ 

(38) a=B^-^B,x\ 

by properly determining the constants A^, A^, B^, B.^. 
To satisfy the differential equations we must have 

while A^ and B^ must satisfy the two equations 

8X^2-1(5X4- 7ya)^2-|w«c = 0, 

- 6X7?, + (5X + 1im)A^ + ^pgac =0, 

which are seen to be identical. The surface conditions give two additional rela- 
tions among the four constants ; these are 

(41) 4X(5„ + i?J _ 2(X + fM)( A> + A) - l(^ + 2^)^, = 0, 

(42) A, + A, + 2B,= 0. 

Solving these equations, we have 

4X + 3/i pga 3X + 2/a pga 

"^iMT+lAp.^jr'^' 2= -2(T9r+ 14m) 7^^' 

7? _ ^^+ ^1^ P9« 7? _ _ 5^+ ^M pga 

""2(19X4-14^) A^ °' '~ 4(19X + 14^) V ^' 

Hence, finally, the values 

(^^) ^ = 2(19xTT47)t'-('^ + ''^)-('^ + '^^*^^T' 

(4*) « = 4(I9X-Vi4/0^'^'^ + ''^^-^'^ + ''^^'^^t'' 

which agree with the known solution for this case. 

(h) Non-gravitating body with given surf ace forces or surface displacements. 
If the two surface stresses are not zero but are known functions of 6 of the 
type above specified, the last members of (32) and (33) become known constants. 
If the surface displacements are given and of the type specified, (32) and (33) 
are replaced by equations obtained in an obvious manner. 

(c) Non-gravitating spherical shell. The solution for a spherical shell may 
be obtained by assuming 

(45) e = A_^x~^-irA_^x-^-irA^ + A^x^, 

(46) a = B_^x-^ + B_^x-^ + B^ + B^x\ 
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the constants being determined to satisfy the conditions at the inner and outer 
boundaries. It may in fact be shown that these equations satisfy (34) and (35) 
if Ag, B^, A^, B^ satisfy (39) and (40), and the remaining coefficients the con- 
ditions 

(47) M_3-(3X + 5M)i?_3 = 0, 

(48) ^_, + 35_, = 0. 

There are thus four constants, available for satisfying the boundary conditions. 
id') hvcompressihle gravitating body. In this case all the gravitational forces 
which are due to changes of density vanish, but the attraction of the surface 
inequality remains. The potential due to this surface layer has already been 
given, being 

This gives in equations (25) and (26) the terms 

Since e^ is a constant, the solution may proceed as in case (a), by assuming 

(49) e = ^0 + ^2^'^ 

(50) a=B, + B^x-. 

In satisfying the equations which determine the constants, special reasoning is 
required by the fact that X is infinite. This presents no serious difficulty, how- 
ever, and will not be here given, since the case of' incompressibility is hereinafter 
treated as a special case of the general problem of a gravitating sphere. To the 
solution of the general problem we now proceed. 

§ 6. Solution of equations in the general case. 

Inspection of the differential equations (25) and (26), and the form of the 
solution in the special cases above considered, suggest the trial of power series 
for the values of e and a . 

The substitution of 

(51) e = ^ cc'% a = B x"^ 

results in the following terms in the two equations : In (25), 

([4(X + ya)m-8M]5„-[t(X + 2M)(m2 + 3m)-4M]^^)x» 
r2(m' + Sm +6) „^ 

(52) +^^4 " 3«K»^?)^(«^™-(- + ^)^™)-^' 



,m + 2 



-l^(2B -A )x'~\. 

5m ^ m m/ J 
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In (26), 

([M(m* + 3m)-2(3X + 5M)]5„ + [(\ + /^)m + (3\ + 5^)]^„)x"' 

(53) +Pff^[- m(J+5) (^^--^'' + ^^'^-^'^'' 

Notice first that the coefficients of x^ in these two expressions are in a con- 
stant ratio independent of m ; hence if m be given a series of values, the coeffi- 
cients of a^ in the two equations may be made to vanish by a single relation 
among the coefficients. 

Again, notice that terms in a;'"+^ result only from the following terms in the 
values of e and a: 

A^ar + A^^,x-+\ B^x"' + B^^,x-+'; 

and that the coefficients of 33™+^ in both equations will vanish if the following 
two equations are satisfied : 

'^^ 3wi(ot-|-5) '- " ^ J n,^ 

lt>.m{m + 7 ) - 6\]5^^, -H [(X -H M)m + 5\ -h 7m]^„^, 

(55) 3 

- pga —. T-^.VoB^ - (to -H 3 )^„ 1 = . 

'^^ m(m-f-6)'- m ^ I y mj 

Finally, notice that the terms of lowest degree in the series for e and a will 
introduce into the two equations terms of like degree which must vanish identi- 
cally. That is, the terms 

Ax"', B^iif 

of lowest degree must be such that the coefficients of cc" in (52) and (53) will 
vanish. The values of m for which this is possible without the vanishing of 
A^^ and B^ are found by equating to zero the determinant of the coefficients of 
A^ and B^ in these expressions ; this determinant reduces to 

|/i(\ + 2/t)(m— 2)m(m -|- 3)(m -|-5). 

The lowest power of cc in the two series must therefore be either 

7?, 03°, x~^ or x~^ . 

Since negative powers of x are excluded by physical considerations (because the 
solution must apply to the point cc = ) , we conclude that the terms available 
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for the solution constitute ascending series of even powers of x beginning with 
the power zero. It is easily found also that we must have * B^=^^Ai^. 
We now assume 

(56) e= A^+ A^x""-^ A^x'+ ■■■, 

(57) a = lA^+B^x'-^ B^x'+.--, 

and determine the relations which must hold among th« coefficients in order that 
the coefficient of every power of x in each of two equations may vanish 
identically. 

From the above discussion it is evident that if »i > , the equations to be sat- 
isfied are (54) and (55), while the vanishing of the coefficient of x^ requires that 



(58) 



[12 " 1 n 

i^o + -5- E - ( 25. - ^.)J = ipgac 

(odd values of m being excluded in the summation). 
The solution of equations (54) and (55) gives 

,rnN A /im^ + 10/im — 6(X — 2/x) n x . -, gap 

^ ' "-+2 /iTO(m + 2)(m + 5)(m + 7) ■- »> V ^ / "'-I X + 2/i 

(60) B,,= S--- o V?-"^ ^^^^%-^ [6^ _ ( m + 3 ) ^ l -^ , 

^ ' '"+2 /im(m + 2)(??i + 5)(»i + 7) L »> v ^ ^ ""-IX + 2/x 

(61) 65 ^,-(m + 5)^ ^,= '^^-^^-^.[QB -{m + 2,)A ^ :c--%- ■ 

\ J m+i \ ^ J m+2 m[m + b)^ " ^ '^ '"-IX + 2/i 

By successive applications of the last equation may be derived the following 
values of the coefficients of x^ and higher powers, n having values 0,1,2, • • ■ : 

(62) A..= -A(65,-5^,)(-^J"[4(n2-f 7.-f9)-6^] ^ . ^ -^l^^^^^ . 



(63)5,,,=2V(65-5^j(-^|^y'|2(rt+3)+(2n+7)^]^ 



n+3 



■9--.(2w-f9)' 
so that the values of e and a may be written 

e=^„ + Ax^-f ^5(65,-5^,)^^^^^ 
(64) 

*This Is evident algebraically. It is also evident otherwise that if the solution has the 
assumed form, B^ must equal J^o ; for Jq and B„ are the values of e and a at the center, where 
these quantities will be related as in the case of a homogeneous strain. 
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a = lA^ + B,x' - ^\ {QB, - 5A,) .-^- X* 



(65) 



X + 2/t 

i:(-ir[2(n+3) + (2r^ + 7)^]^ 



w + 3 

9 ...(2« + 9] 



( i^^A" 



in which A^, A^, B^ must satisfy (58) and the surface conditions (32) and (33). 
For substitution in (58) we find (for n = 0, 1, 2, • • •) 



(66) 



4+ 2re 






71 + 3 



9 ••• (2W + 9) 



4(.^+7. + 9)-6^]^-^^2^^^^., 

w + 3 



2(?i + 3) + (2w + 



In order to express as concisely as possible the values of e and a and the 
three equations to determine the constants, the following notation may be 
employed : 

(67) X(.) = t(-l)"r-- . " . "^ -^^1 '^ + ^ 



(68) M{z) = j:^{-1Y 



(69) N{z) = ±{-\f 



(70) L\z) = ±{-ir 



(71) M'{z) = ±{-iy 



W + 6 + 



\ 1 »i + 3 

^j7-9 ... {iiV^rrf^' 



9 ... (2w + 9) ' 



4(n2 + 7»+ 9 



)-^]l 



X-| (n + 2)(7i + 3) ^ 



9 ...(2« + 9) ' 
i(,^ + 3) + (2n + 7)-J^ g __^2n^gy , 



(72) 



6 = 



gap 



X + 2/t' 



The solution then takes the following form, the last two of the five equations 
being derived from the surface conditions (32) and (33) : 



{A) 
(C) 



4(-lo+A) + 



e = A^ + A^x? + 2^(6^2 - ^A)hx*L{bx^), 
a = i^„ + B.x' - 2\{6B, - 5A^)hx*M{bx'), 
140 



A 



Sb{X/fi+2] 
+ (6Z^,-5A)(2-3^^--^-62N.(6))+-V-c = 0, 
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[2 + r ^ + 2) sl {A^ + A,) + \A, - (6^, - hA,) 

(E) 3(^. + 4,) + 5^, + (6B,-6^)[l-g[i(6)-2Jf(6)]]-0. 

The functions represented by the infinite series may be otherwise expressed, 
depending ultimately upon the exponential integral 






^^'^'^dx . 



For the purposes of computation, however, it is doubtful whether they can be 
reduced to expressions more convenient than the series, which converge with 
satisfactory rapidity unless 6 is relatively large. In applying the results to a 
body comparable in size and mass with the earth, for example, it can hardly be 
supposed that h is as great as 2 ; and for 6 = 2 it is found that eight or ten 
terms of the series suffice, even in case of the function X', which is the least 
rapidly convergent of all. 

§ 6. Modified procedure for solving equations. 

It may be well to point out that by changing the procedure the algebraic 
work involved in the solution of equations (25) and (26) can be somewhat 
abridged. 

The equations of equilibrium of a volume element may be abbreviated by 
introducing the elementary rotations t<r^, tjr^, w^. The general forms of these 
equations are given in Part IV of this paper, equations (117), • • • , (123). "When 
the strain is symmetrical about the initial line of the system of coordinates, 
vy^ and m-g vanish, and the equations reduce to the two following, which are iden- 
tical with (1) and (2); 

aA 2a 5 
(X + 2,.)^-^-:^^^(^,sin^) + pi^,= 0, 

1 5A 2m 5 

(X + 2,.)^^ + ^-^.^^-^-C'-*sin^)+pF,= 0. 



If we now write 
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and introduce everywhere y and z instead of e and a, equations (25) and (26) 
become 



(\+2m)2/ + m 



dx 



"''^"[5^ J y'«''f^a'' + -^j ^ + (|ei+ic)x2j=0. 



It is easily shown that these are satisfied by series of the following forms : 

y = C^x^ + C.jc* + : • • , z= D^x? + Z»,x* + • • • , 
in which the coefficients satisfy the equations 

C„ m(w + 5)' X + 2yLt' 



Z> 



/'5'« 



4a — 2e — |x T- = ^ 
<* ax 



(7„ »z(m + 5) M ' 

and an additional relation obtained by equating to the coefficients of x? after 
substituting the values of y and z in the differential equations. One coefficient 
thus remains arbitrary.* 

It is now possible to determine e and a from the equations 

-y, f(e-2a-x£)=.. 

Thus if we assume 

e= 2^ x"% a = ^B x"", 
we find 

(m + 3)J,„-65„=-|C„, 

^„-(m + 2)5„ = |Z»^, 

from which A^ and 5„ are determined in terms of C and D (and therefore 
in terms of the one arbitrary constant above mentioned), except in the case f 
?w = . In this case we may take J.^ arbitrarily and make 



* Another solution consists of ascending series beginning with C—ss;""' and Z)— 3^~' ; but nega- 
tive powers of X are not permissible in the physical problem. 
t The case to ^ — 5 is excluded by physical considerations. 
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That is, the values of e and a contain the same powers of x as the values of y 
and 2, together with the absolute terms A^ and J^q. Two arbitrary constants 
remain for satisfying the surface conditions. 

Although this process results in some saving of algebraic work, it does not 
appear to suggest any shortening of the labor involved in numerical computations 
such as those given below. 

It is of some interest to note the form of the differential equation resulting 
from the elimination of z : 



^^y\ o^'^y a,. , w« 



+ 2a'-/.-62/ + 



d7^^ dx ^^X+2fi 



(.3| + 6.V) = 0. 



We obtain also the equation 



, dh „ dz 



^dz ^ dz „ pqa „ 

x^ -j-„ + 2x , 6g — '-^- a?y = ; 

ax^ dx fi -' 

but the complete elimination of y obviously leads to an equation of the fourth 
order in z. 

§ 7. External potential and inertia-constant. 

One of the quantities which indicates the yielding of the body is the change 
in the external potential, or (what amounts to the same thing) in the difference 
between the principal moments of inertia. 

Let / denote the moment of inertia about a diameter before strain, O and A 
the principal moments after strain (the former being taken about the axis of 
symmetry), and let 

(73) /=^. 

The value contributed to C — A by an elementary shell of thickness dr and 
density 

- pA = - p(^4a - 2e - fr ^^^P, 



IS 



l-Trr^' 



V(^4«-2e-trJ)c?r, 
while the part contributed by the surface inequality is 

hence 

(74) C-A=^^ j^(^4.a-2e- |r J) r^dr + J^-rrpa^e,. 

Integrating by parts the term containing de/dr, we obtain 
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(75) C-A = ^~^ r{Sa+e)r*dr=^^^ C {Sa+e)x'dx, 
and 

(76) /= ^~^ = 2 r'(3a + e)xMx. 

This accords with the result already given for the potential due to the change 
of configuration. Thus the external potential is known to have the value 

y{C-A) 

while it is also represented by the first term of the value of V given in (20) if 
the upper limit of integration be a ; i. e., we must have 

yjC-A) 4gP, r" „ IQ-TrpyP, r 
p ^2=-'^^J {^a+e)r dr = i[^^I~ J {^cc + e)r'dr, 

which agrees with the result given above. 

Using the values above found for e and a, we may write 

('^7)/=^~ = A + A + (6i?.-5A)[7 + 4(2G^(&)-3^i^(i))], 
in which 

(78) •p(')-l:(-ir,-T9T^f2-'-^j-. 

§8. Case of incompressibility. 
In the case of incompressibility X ^ oo , 6 = 0. Putting b' = gap/ fj. , we find 

(80) »^(»)-x?^-,(^')-0-A»-. 
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Equation ( C) now requires that 

(81) 65,- 5^,= 0, 

(82) (65,-5^,)~=|'[4^„ + 4(l+g)^,+ V^c]- 
while (Z>) and {E) give 

(83) {QB,-5A,)\ = {2 + l')A^ + {& + V)A„ 

r 

(84) ^ 3J„ + 8^, = 0. 
Solving, we obtain 

.„,, J _ 8^' , _ 35' _ 55' 

(^80) ^0-38 + 45''^' ^2--38^4^,'C, ^2 " ~ 2(38 + 45' ^' 

while all other coefficients in the values of e and a vanish. We thus have 

(86) h38T45'(«-^-^)' 

(^^) ^ = 2(38 + 45') (^-^-^)- 

which agree with the known solution for the case of incompressibility. 
We find also 

(88) /=A + A = 3-^" 

which is identical with the surface value of e, as should be the case when 
\ = 00. 

§ 9. Case of zero rigidity. 

This is the only case in which the solution fails. If /^ = , equations (25) and 
(26) cannot be satisfied unless A = ; and if A = 0, the surface condition (32) 
cannot be satisfied unless X = 00. 

§ 10. Numerical results. 

In applying the formulas above found, it is convenient to choose a definite 
value of \//A and make computations for a series of values of 5 or gpa/( X + 2/u.) . 
Some of the results thus found are given in the accompanying tables and dia- 
grams,* which for the most part are self-explanatory. The quantities of par- 
ticular interest are the surface values of e and a and J". As an aid to the 



'^In all the diagrams it is to be understood that ordinates are multiplied by c. 
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verification of the computations, the values of the functions L, M, etc., are 
included in the tables. 



Table I. Results for the case Xj ii = 1 . 

(All numbers below the horizontal line are to be multiplied by c. Numbers in square 
brackets [ ] obtained by interpolation ; all others computed. ) 



. gap 


2 


1.8 


1.6 


1.4 


1.2 


1 


.8 


.667 


ft/gap 


.167 


.185 


.208 


.238 


.278 


.333 


.417 


.500 


Lib) 


.8946 


.9347 


.9773 


1.0224 


1.0703 


1.1212 


1.1753 


1.2133 


J>Hb) 


.4618 


.4746 


.4879 


,5018 


.5163 


.5315 


.5474 


.5585 


N(b) 


.2579 


.2641 


.2706 


.2774 


.2844 


.2917 


.2994 


.3026 


L'{b) 


1.4000 


1.4978 


1.6042 


1.7194 


1.8444 


1.9801 


2.1273 


2.2326 


M'{b) 


.7987 


.8319 


.8668 


.9039 


.9434 


.9853 


1.0298 


1.0610 


F{b) 


.03798 


.03880 


.03964 


.04052 


.04142 


.04236 


.04333 


.04400 


G(b) 


.02999 


[.031] 


[.033] 


[.034] 


[.036] 


.03758 


.03936 


.04061 


6Bj— 5^, 


3.710 


2.857 


2.141 


1.552 


1.067 


.7004 


.4152 


.2701 


^o=''o 


.5146 


.5251 


.5218 


.5058 


.4758 


.4327 


.3760 


.3303 


A, 


—.5752 


—.4954 


—.4228 


—.3570 


—.2964 


—.2405 


—.1883 


—.1551 


B, 


.139 


.0633 


.0045 


—.0388 


—.0677 


—.0837 


—.0877 


—.0842 


«i 


.216 


.230 


.238 


.241 


.237 


.225 


.204 


.1843 


/ 


.436 


.415 


.390 


.361 


.327 


.289 


.246 


.213 


«! 


.253 


.224 


.196 


.169 


.142 


.117 


.0927 


.0767 



Table II. Results for the case Xj ii — 2 . 

(Numbers below the horizontal line to be multiplied by e. Numbers in square brackets 
obtained by interpolation ; all others computed. ) 



X+2|it 


2 


1.8 


1.6 


1.4 


1.2 


1 


.8 


.667 


/i/gap 


.125 


.139 


.156 


.179 


.208 


.250 


.3125 


.375 


L{b) 


.6667 


.7019 


.7394 


.7793 


.8218 


.8671 


.9153 


.9493 


M{b) 


.7117 


.7312 


.7516 


.7729 


.7952 


.8185 


.8428 


.8597 


N(b) 


.2959 


.3029 


.3103 


.3179 


.3258 


.3341 


.3427 


.3467 


L'{b) 


.9922 


1.0771 


1.1697 


1.2706 


1.3806 


1.5005 


1.6311 


1.7367 


M'(b) 


L2321 


1.2829 


1.3365 


1.3933 


1.4537 


1.5179 


1.5860 


1.6328 


F{b) 


.03798 


[.0388] 


[.0396] 


[.0405] 


[.0414] 


.04236 


.04333 


.04400 


0(6) 


.02999 


[-031] 


[.033] 


[.034] 


[.036] 


.03758 


.03936 


.04061 


6.Bj — 5^2 


3.749 


2.915 


2.209 


1.623 


1.142 


.7577 


.4596 


.3056 


■^o = f<, 


.5537 


.5852 


.5984 


.5935 


.5698 


.528 


.4666 


.4144 


A, 


—.4655 


—.4309 


—.3938 


—.3542 


—.3118 


—.267 


—.2195 


—.1862 


s. 


.2369 


.1267 


.0400 


—.0247 


—.0695 


—.0962 


—.1063 


—.1042 


«i 


.2964 


.3078 


.3135 


.3131 


.3049 


.289 


.2611 


.236 


/ 


.519 


.496 


.469 


.438 


.401 


.358 


.307 


.269 


«! 


.2914 


.259 


.228 


.199 


.170 


.142 


.1140 


.096 
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Table III. Results for the case X//ji, = 3 . 

(Numbers below the horizontal line to be multiplied by e. Numbers in square 
brackets obtained by interpolation ; all others computed. ) 



223 



A +2(1 


2 


1.8 


1.6 


1.4 


1.2 


1 


.8 


.667 


fi/gap 


.10 


.111 


.125 


.143 


.167 


.200 


.250 


.300 


L(b) 


.4388 


.4691 


.5016 


.5362 


.5733 


.6129 


.6553 


.6853 


M{I>) 


.9616 


.9880 


1.0154 


1.0441 


1.0740 


1.1054 


1.1381 


1.1610 


N{b) 


.3329 


.3417 


.3499 


.3584 


.3672 


.3764 


.3861 


.3906 


L'{h) 


.5840 


.6566 


.7354 


.8218 


.9168 


1.0212 


1.1349 


1.2414 


M'(b) 


1.6654 


1.7339 


1.8060 


1.8826 


1.9641 


2.0505 


2.1423 


2.2046 


F{b) 


.03798 


.03880 


.03964 


.04052 


.04142 


.04236 


.04333 


.04400 


G{b) 


.02999 


[.031] 


[.033] 


[.034] 


[.036] 


.03758 


.03936 


.04061 


6^2 — 5^2 


3.910 


3.032 


2.320 


1.721 


1.224 


.823 


.497 


.343 


^0 = ''o 


.545 


.598 


.643 


.654 


.642 


.605 


.532 


.489 


A, 


-.340 


— .348 


— .351 


— .339 


—.317 


— .285 


-.239 


— .214 


B, 


.368 


.215 


.094 


.0043 


— .0602 


— .1003 


-.1163 


— .1215 


fi 


.3487 


[.365] 


.370 


.369 


.360 


.341 


[.314] 


.281 


f 


.580 


.556 


.533 


.501 


.463 


.417 


[.370] 


.318 


«i 


.327 


.290 


.258 


-.226 


.195 


.164 


.131 


.112 



In Figs. 2 and 3 are shown curves representing e^, a^ and f as functions of 
fi/gpa, for certain constant values of X//K.. It is seen that, as X/yx increases, 
the curves for e, and f both approach the curve (A) which represents the case 
X = 00. For a given value of fj, both gj and _/ increase with X ; in other words, 
the effect of compressibility is to decrease both the ellipticity and the external 
potential corresponding to a certain rigidity. This effect is much more impor- 
tant upon e, than upony. Thus if /i = .S5pga, a variation of X//u. from 1 to 
oo would change f/c only from .279 to .289, while e^/c would range from .221 
to .289. For smaller values of /i the effect of compressibility is greater. An 
interesting feature is the fact that for a given value of X//i there is a maximum 
value of Cj ; thus if X=/i, a decrease of /jl below about .24pga causes e^ to decrease. 

Fig. 3 shows also that the value of otj is increased to an important extent by 
compressibility. 

II. Strain of an elastic sphere covered by a shallow ocean. 
§ 1. Statement of problem. 

In applying the theory of the deformation of a sphere to the case of the earth, 
an important question is the effect of the ocean. Some light on the order of 
magnitude of this effect is thrown by the solution of the following problem : 

To determine the strain of a gravitating elastic sphere, the surface of which 
is covered by a shallow ocean in equilibrium. 

The solution above given may without difficulty be extended to this case, as 
follows. 
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§ 2. Forces due to ocean inequality. 

Let c'p = density of ocean, 

e = ellipticity of ocean surface. 

Then there must be introduced into the solution the forces due to an ocean 
layer of thickness 



-|(e-eJaP, 



.60 



.50 



.40 



.30 



.20 



.10 





\a 














^ 










/■ 


^?^ 


^ 










^^5^ 


::::::^ 


^**^ 




V-"=J 



























.10 -30 ij./gpa -80 ' ' .40 .50 .60 

Fig. 2. Curves showing efiEect of compressibility upon surface ellipticity. 

These forces wUl be (a) the bodily force due to the attraction of the layer, and 
(b) the surface pressure due to the weight of the layer, 
(a) The attraction is given by the potential 



(89) 



- = - ^7r7C>(e - e,yP, = _ gc'(e - eJr^P,, 
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which adds to (23) and (24) the terms 

d V 4:QP 



r 06 oa ^ ^' 
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Fig. 3. Curves showing effect of compressibility upon angular displacement 
and external potential. 

and therefore to (25) and (26) the terms, 

(90) _4^pac'(e-ejx% 

(91) ^gpac'{e—e^)x\ 
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(b) The weight of the ocean inequality per unit area is 

hence there must be added to the second member of (32) the term 

(92) ygac'{e-e,). 

§ 3. Value of elUpticity of ocean surface. 

The value of e (the ellipticity of a surface which is level under the action of 
gravity and the disturbing forces) may be expressed in terms of f, e^ and c , as 
follows : 

For the solid sphere the value of C — A is If, while for the ocean inequal- 
ity it is 7c'(e — e^). Hence by a familiar result in the theory of potential, 

(93) I[f+c'{e-e,)]=iMa\e-ic). 
Or, since for a sphere of uniform density /= |Jfa^, we have 

(94) /+c'(e-eO = f(e-ic), 
whence 

^^^-' ^- 10 -6c' ' 

and 

The terms to be added to (25), (26) and (32) thus become 

(98) ,,„[... ^iL+^f..]. 

r ,6/-10e, +5c'l 

§ 4. Statement of solution. 

The values of e and a for this case are expressed by the same equations [A) 
and (5) as before, but two of the three equations (C), (D) and (^) which 
determine the constants must be changed as follows. 

In (C) replace -^-c by [50c + 6c'(6/- 10eJ]/[3(5 - 3c')]. 

Add to the first member of (D): 



(^^) 



c'(5c + 6/"-10ej) 
10^^67 ' 
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Equation {E) remains unchanged. 

Equations ( (7) and ( Z) ) are thus replaced by the following : 

140 

12c' 20c' 50 _ 

+ 5^^Ji'f- S^^'^i + 3(6-3c')''= ' 

[2 + (^+2)5]^,+ [6 + ^ + 2)5] ^,-(65,- 5^,) 

(A) x{^-^,[(3^+2+Q+2)6)x(6) + 6^i.f(6) + 2Q+2)x'(6)]} 
/ ^ ^\ , / 3c' ^ 5c' 5c' \ „ 

When f and e^ are replaced by their values 

(100) e^ = A, + A,+ Jj(65, - bA,)hL{h), 

(101) /=^„ + ^, + (65,- 5^j[^+ 5^(^2(^(6) -3^i^(6))], 
the equations become 

20(1 — c') 140 

-T3^ (^0 + ^^) + 3T(x7^+2) ^1^ + ^^^^ - ^^^) 

(^') X [2-3T^2l-^^(^)+54W(Hi^G^(^)-|>^(^)-/.*^(^))] 

50 

+ 3(6-3c')''~ ' 

[2 +(^+ 2)5-^1^3^9] (^„ + ^,) +4^,- (6iS,- 5^,) 

>< j^-M[(3^2 + (^+2)6)5Z(6)+6^6Jf(6)+2(^+2)6i'(6)] 

/ >- ^\ , 5c' 
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To these we add equation (E) unchanged : 

(B') S{A,+ A,) + 5A, + {QB,-5A,)^l + l{L{h)-2M'{h))j = Q. 

These three equations determine the constants A^, A^, B^. 

§ 5. Value of inertia-constant for comhined sphere and ocean. 
Let/" represent the value of (C— A) /I for the whole system, then 

/=/+c'(e-e,), 
which by (96) reduces to 

(102) / ^s'-Tg^'/- 5^37^- + 2(5 -3c') ^- 

§ 6. Case of incompressibility. 

If X = 00, equations (^A) and {B) hold, and equations (C) and (Z>') require 
that 

(103) _ 6B,-5^, = 0. 
Also, on writing 

(104) \{%B,-bA,)^C, 
(C), (Z)') and {E') become 

(105) C-7A-'|~^&'(A + A)-2(5^«=^' 

(106) C-4^,-(2+^|if^U')(^„+^,)-2-(-|^^c = 0, 

(107) 3(^„ + ^,) + 5^, = 0. 

The solution of these equations gives 

20(1 -c')5' 

» ~ 19(5^3 3c'y+T0Xl -c)h''^' 

n08^ A- - _i5(l-0i^ , 

'^ -' ^~ 38(6 - 3c') + 20(1- c')6'' 

p 25(1 -c')6' 

^ ~ 76 (5 -3c') 4-40(1 - c')!)'^' 

If c' = , these reduce to (85), as of course they should. 
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§ 7. Numerical results. 

For the purpose of showing the order of importance of the ocean effect, com- 
putations have been made with a value of c equal to the ratio of the density of 
sea-water to the mean density of the earth. This value is taken as 5/27. 

The computations made cover the two cases \//a = l, \/fj, = 00, and the 
results are shown in the accompanying table and diagrams. 

For the case of incompressibility 

226' 22 

114 + 116 114JL+11 

gpa 

(110) '- = f ^(1-K), 

114^- + 11 
gpa 

(111) " = ^ (1-K)> 

114-^+11 

ffpa 

(112) J = 7 = '(^o + A) = 4sr^- 

11 gpa 

The last expression is written in a form for comparison with the value without 
ocean effect: 

e, f 5 



(113) 



' 38^ + 4 



gpa 
From (102), 

(114) /'=ii«i + lV- 

Values of e^, f and _f for this case, also corresponding values without ocean 
effect, are given in Table IV, while the same results are shown graphically in 
Fig. 4. 

For X = /a , computations have been made by means of equations {A'), {B'), 
{ C), (D'), (^'), with the results given in Table V and shown graphically in Fig. 5. 

Inspection of the curves shows that the effect of the ocean on the strain of 
the sphere is quite appreciable even in the case of incompressibility, and is very 
considerable in the case \ = /a . A fact of equal importance in the interpreta- 
tion of observational results as to the yielding of the earth is shown by the 
curves representing y. These are nearly identical for the two cases X = 00 , 
X = yu. (Figs. 4 and 5), and in both cases are very considerably above the curves 
representing f. It thus appears that an important part of the change of 
external potential is due to the ocean — a result which will be referred to again 
in the following discussion of the yielding of the earth. 
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Fig. 4. Curves showing efleot of ocean in case /I = oo. (Fq and F^ with ocean, F without. ) 
Table IV. Effect of ocean in case X= oo. 





Without ocean. 


With ocean. 
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oi 
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c 
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.375 


1.250 


.375 


1.250 


1.250 


.05 


.254 


.848 


.247 


.823 


.858 


.10 


.193 


.642 


.184 


.613 


.666 


.15 


.155 


.516 


.146 


.488 


.551 


.20 


.129 


.431 


.122 


.407 


.477 


.25 


.111 


.370 


.104 


.348 


.423 


.30 


.097 


.325 


.091 


.304 


.383 


.35 


.087 


.289 


.081 


.270 


.352 
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.078 


.260 


.073 


.243 
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.071 
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.066 
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Fig. 5. Cnrves showing effect of ocean in case X^ii. (Curves E, F without ocean ; E^, i^», 

i^Q with ocean.) 



Table V. Effect of ocean in case \ = ;tt . 







Without ocean. 


"With ocean. 
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»i 


ei 


/ 


»! 


ei 
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/' 
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c 


c 


c 


c 


c 
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2 


.167 


.253 


.216 


.436 


.257 


.119 


.379 


.505 


1.8 


.185 


.224 


.230 


.415 


.229 


.145 


.365 


.485 


1.6 


.208 


.196 


.238 


.390 


.198 


.163 


.345 


.458 


L4 


.238 


.169 


.241 


.361 


.170 


.176 


.322 


.429 


1.2 


.278 


.142 


.237 


.327 


.143 


.181 


.294 


.397 


1.0 


.333 


.117 


.225 


.289 


.118 


.179 


.263 


.363 


.8 


.417 


.093 


.204 


.246 


.094 


.167 


.225 


.322 


.667 


.500 


.077 


.184 


.213 


.077 


.151 


.194 


.291 
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III. Effect of compeessibility upon estimates of the 

RIGIDITY OF THE EARTH. 

The estimates of the rigidity of the earth, made originally by Lord Kelvin 
and since by others, are based upon a comparison of what is known regarding 
the actual yielding of the earth under known disturbing forces with the computed 
yielding of an ideal body under like forces. This ideal body is a homogeneous, 
elastic, incompressible sphere of mass and size equal to those of the earth. It 
is, of course, not to be expected that such a comparison can lead to anything 
more than a rough average value of the modulus of rigidity of the material of 
the earth, since the actual earth is certainly far from uniform in density, and 
probably far from uniform in elastic properties. It is, however, a matter of some 
interest to notice in what way these estimates are modified when compressibility 
is taken into account. In answering this question it is needful to refer briefly 
to the present status of the evidence regarding the actual yielding of the earth.* 

§1. Nature of evidence of yielding of the earth. 

Evidence regarding the actual yielding of the earth relates to the two quanti- 
ties above represented by e^ and/. The former is the surface ellipticity due to 
strain, the latter measures the change in the gravitational potential at the earth's 
surface due to the strain. The number f is, in fact, the value of (C — A) /I 
due to the strain. It may be remarked that for a homogeneous incompressible 
body e^ and f are equal, but for a compressible body f depends upon the values 
of both e and a throughout the body. 

There are three lines of evidence of yielding of the earth to tidal or centrifugal 
forces : (a) Tidal observations, (6) experiments with the horizontal pendulum, (c) 
the movement of the axis of rotation within the earth as determined by latitude 
observations. 

(a) Nature of evidence from tidal observations. The height of an equilib- 
rium-tide at any locality can be computed on the assumption that the solid earth 
is unyielding, also on the assumption thaty and e^ have any definite values. A 
comparison of these results with the observed tidal oscillation gives eviderice as 
to the actual values of f and e^ . It may, in fact, be shown without difficulty 
that if the attraction of the ocean itself be neglected, the ratio of the actual 
fluctuation of the water relative to the solid ground to the fluctuation on the 
assumption of no yielding of the solid earth is 



* A clear and non-technical summary of this evidence is given by G. H. Daewin in a paper 
entitled " The Rigidity oftlie Earth," in Eivista di Scienza, vol. 5 (1909). See also a paper 
by A. E. H. Love entitled " The Yielding of the Earth to Disturbing Forces," in the Proceed- 
ings of the Eoyal Society, Series A, Vol. 84. 
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c c c 

It is here assumed that the tide has its equilibrium value, an assumption which can 
reasonably be made only in the case of tides of long period. It is therefore only 
the fortnightly and monthly tides that have been available for this comparison. 

It should be said also that for full generality the term fjo must be multiplied 
by HjMa}. For a homogeneous sphere the value of this is 6/5, but for the 
earth it is very nearly 1. It is the actual value for the earth that must be used 
in the application made below in § 2. 

(5) Nature of the evidence from experiments with the horizontal pendulum. 
If the change in the direction of apparent gravity relative to the earth at any 
locality due to tidal forces can be measured, a comparison of this effect with the 
effect computed on the assumption of absolute rigidity of the earth will give 
evidence as to the actual yielding of the earth. The actual effect is so minute 
that it must be greatly magnified in order to become measurable. This magni- 
fication has been accomplished by means of the horizontal pendulum. It is not 
difficult to show that the actual deflection of such a pendulum due to tidal forces 
bears the same ratio to the deflection that would occur if the earth did not yield 
as the amplitude of an equilibrium tide bears to the amplitude it would have 
if the earth were unyielding. The quantity determined from horizontal pendu- 
lum observations is therefore identical with that inferred from tidal observations. 

(c) Nature of the evidence from the periodical variation of latitudes. The 
periodic motion of the axis of rotation within the earth furnishes evidence as to 
the yielding of the earth to centrifugal forces. This evidence is of a different 
character from that obtained from the tides or the horizontal pendulum, in that 
it leads to an estimate of the value of fjc instead of to a relation between fjc 
and Cj/c. The change in the value of (O— -A)/! due to the yielding of the 
earth to given centrifugal forces may be inferred directly from a comparison of 
the actual period of the free nutation of the earth's axis with the period com- 
puted for an unyielding body having the earth's actual figure. It may, in fact, 
be shown that the actual nutation period is to a close approximation the same 
as that of an unyielding body having the figure which the earth would assume 
if centrifugal forces were annulled.* 



*This was shown by S. S. Hough (Philosophical Transactions of the Royal 
Society, ser. A, 1896) for the case of a homogeneous, incompressible elastic spheroid, and has 
been proved by J. Larmor and others to hold without other restriction as to the character of the 
body than that its density and elastic properties have approximate spherical symmetry. A proof 
of this was given by the present writer in a paper presented to the American Mathematical 
Society, San Francisco Section, Dec. 20, 1902, of which only an abstract was published, 
Bulletin of this Society, Vol. 9 (1903), p. 299. The reasoning employed was identical with 
that published by Larmor, Proceedings of the Royal Society, ser. A, vol. 82 (1909). 

Trans. Am. Math. Soc. 16 
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§ 2. Numerical results expressing actual yielding of the earth. 

The actual value of (C— A) 1 1 for the earth is 1/308, so that the period of 
free nutation if it were absolutely rigid would be 306 sidereal days. The period 
of the movement of the axis about its mean position, as inferred from the obser- 
vations of the International Latitude Service carried on systematically since 
1899, has varied from 427 days to 442 days. The former gives 

/=a6-dT = -00093; 
the latter, 

/=^i6-?i2 =-00100. 

Since the ratio of equatorial centrifugal force to gravity is c = 1/289, the cor- 
responding values of fjc are .269 and .289. The former of these values has 
hitherto usually been accepted. It is nearly equal to 4/15. 

The first analysis of tidal observations, with reference to the question of the 
yielding of the earth, was made by G. H. Darwin, the data and results being 
given in Thomson and Tait's Natural Philosophy. From all the data then 
available which were regarded as trustworthy, the value found for the ratio of the 
actual tidal amplitude to that computed on the assumption of an absolutely rigid 
earth was .676 ± .076 ; but a separate analysis of the more consistent observa- 
tions gave .931 =h .056. The result was not regarded by the author as conclu- 
sive.* Analyses since made by others have been confirmatory of the former of 
these values, and about the same result is inferred from observations with the 
horizontal pendulum, f 

If this result be accepted, we have (very nearly) 

(115) 1+Z_2'' = ^. 

^ ^ ceo 

And if in this we substitute 

c ~ 15' 
there results 

e,_ 3 
c ~ 10 ■ 

§ 3. Comparison of actual yielding with that of an ideal elastic solid. 

It is pointed out by LovE^ that the above values oi f/c and e,/c cannot 
coexist in the case of a homogeneous, incompressible elastic solid having the 
earth's size and mass, and it is easily seen from the results found in the first 

* "We thus see that the more consistent observations seem to bring out the tides more nearly 
to their theoretical equilibrium values with no elastic yielding of the solid." 6. H. Dae win, 
in Thomson and Tait's Natural Philosophy, Part II, p. 460 (edition of 1890). 

tW. SCHWEYDAE, Beitrage zur Geophysik, vol. 9 (1908) ; O. Hecker, Veroffentlichung 
des Kgl. Preussischen Qeodaiischen laslituies, No. 32, Berlin (1907). 

tLoo. cit. 
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part of this paper that the discrepancy is increased rather than diminished if 
the effect of compressibility is considered. This is shown clearly by Fig. 6, in 
which the curves represent the relation between e^jc and fjc for several values 
of X/jtt, while the straight line AB represents equation (115), or 
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Fig. 6. Curves showing relation between p^ and / for several values of ? //i . 
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As X//i increases, the curves approach the limiting line e^ =-f\ and until X//it 
reaches a value considerably greater than 3 , the curves will fail to intersect AB . 
It is obvious that if X/jit be taken great enough so that (116) can be satisfied, 
the value oifjc will be greater than ^^^ ; being in fact 3 for X = 00 , and greater 
than J for any finite value of X//u,. 

It may be remarked that the discrepancy between' 
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is after all not remarkably great, considering the great delicacy of the observa- 
tions upon which equation (116) is based. It may perhaps be still open to 
question whether the accepted values are so thoroughly reliable as to exclude the 
possibility of errors of the magnitude of these discrepancies.* 

It may be noted, also, that a value of 1 +fjc — 2ej/c much less than the 
greater of Darwin's results would remove the discrepancy. Thus the following 
values would be consistent for the ideal homogeneous sphere : 



A 
I* 


t 
c 


c 


e c 


00 


.270 


.270 


.730 


3 


.270 


.245 


.780 


2 


.270 


.235 


.800 


1 


.270 


.220 


.830 



Although the value of the surface ellipticity accepted as most reliable is more 
nearly consistent with the computed value for an incompressible than for a compres- 
sible sphere, it can hardly be supposed that this furnishes any important evidence 
as to the actual compressibility of the materials composing the earth. As stated 
above, only a roughly approximate agreement is to be expected between the behavior 
of the actual earth and that of an ideal body assumed to be uniform in density and 
elastic properties. Moreover, there are certain questions not generally considered 
wnicn appear to require attention in any complete study of the theoretical basis 
of the comparison. One of these questions relates to the effect of the ocean. 

* The most elaborate and apparently the most trustworthy series of observations heretofore 
made with the horizontal pendulum are those of Heckeb (loc. cit.). A remarkable feature of 
the results obtained for the lunar effect upon the yielding of the earth is the variation of this 
effect with the azimuth of the moon. The value found for the ratio of the actual deflection of 
gravity to the deflection computed for a rigid earth is a minimum when the moon is on the 
meridian, and increases continually with the azimuth. The minimum value of the ratio is about 
0.36, for azimuth 45° it is about 0.50, and for azimuth 67° about 0.90; and although as the 
azimuth approaches 90° the ratio becomes indeterminate, it appears to approach unity. No 
satisfactory explanation of this variation seems to have been suggested. If it is accepted as 
proving that the yielding of the earth varies greatly with the azimuth of the disturbing body, it 
appears to render of little value conclusions based upon the ordinary theory which assumes that 
the yielding is symmetrical about the axis of symmetry of the disturbing forces. At all events 
these results can hardly be regarded as establishing any definite value of the ratio of yielding. 
[The above numerical values of the ratio for different azimuths were obtained by measurement 
from the diagram given by Heckkr showing the observed lunar effect and the effect computed 
for a rigid earth ; assuming that the deflection of apparent gravity always agrees in azimuth with 
the moon.] Hsckkb's results are discussed by G. H. Daewin in the paper already cited. This 
paper was received too late to permit here any quotation of suggestive comments. 
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§ 4. Effect of the ocean on estimates of the yielding of the earth. 

The question as to the effect of the ocean involves two parts, — the first, what 
effect has the ocean on the actual yielding of the solid earth, and the second, 
what direct effect has the ocean on the observational results? 

While no exact answer can be given to either of these questions, something 
can be inferred as to the degree of importance of the ocean effect from the 
results reached in Part II of this paper, and shown in Tables IV and V, and in 
Figs. 4 and 5. These results apply to the two cases X = oo and Xf/x = 1 . It 
appears that the effect of the ocean upon the strain of the sphere is quite appreci- 
able in both cases, and in the case X = /u. becomes of considerable importance. 
Although the angular displacement a^ is affected but slightly, both e, and f are 
very considerably changed by the influence of the ocean. 

Moreover, while these quantities are both materially decreased, the total exter- 
nal potential due to the solid sphere and ocean is considerably greater than the 
potential when the ocean is absent. This is true both for X = jj, and for X = oo , 
as is seen from the values oi f'/c. (It is of interest to note that the curves 
showing the relation between f'fc and i^'/gpa are nearly identical for the two 
cases.) 

In estimating the value of yu. from the value of (C — ^)/Z inferred from the 
prolongation of the period of free nutation, it is therefore important to know 
whether the quantity determined applies to the solid earth alone or to the com- 
bined solid earth and ocean. Thus if it applies to the solid earth alone, the 
value f/c = 0.27 would correspond to 

— = about 0.37, 
gpa 



while if 0.27 is the value of f'/c, we find 



— = about 0.57. 
gpa 

The former would give /m = 1.30 x 10'^ while the latter would give /u,= 2x 10". 
These values are in C. G. S. units (dynes per square centimeter), the value of 
gpa in these units being taken as 3.5 x 10'^. 

Whether the change in the value ot (C — A )/I due to centrifugal forces, as 
inferred from the nutation period, applies to the solid earth alone, or whether 
an important part of it is due to the change in configuration of the ocean, is a 
question which must be answered in the light of the theory of the effect of 
yielding upon the nutation period. It seems a reasonable assumption that, for 
a disturbance for so long a period as 427 days, the ocean always has its equi- 
librium figure. If this be assumed in applying the equation of angular momen- 
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turn,* it is difficult to avoid the conclusion that the change in the value of 
(C— -^)/I inferred from the prolongation of the nutation period includes the 
yielding of the ocean; in other words, that the number 4/15 is the value of 
j'/c rather than oifjc. 

What effect the ocean may have upon observations with the horizontal pendu- 
lum is a more difficult question, because the equilibrium theory must be dis- 
carded in dealing with tides of short period. 

IV. Strain of a compressible, gravitating elastic sphere due to 

DISTURBING FORCES OF WHICH THE POTENTIAL IS A SPHERICAL 
HARMONIC OF ANY DEGREE NOT LESS THAN THE SECOND. 

The following solution of this problem follows the general lines of the solu- 
tion already given for the case in which the degree of the harmonic is 2. It in 
fact includes that as a special case. In comparing the algebraic results it must 
be noted that there is a slight change in the notation ; the quantities e and a as 
used below must be multiplied by the factor — | in order to agi'ee with e and a 
in the solution given above. 

§ 1. Physical analysis of the prohlem. 

The physical analysis is identical with that already given in the special case, 
and need not be repeated. 

§ 2. Formation of the equations of equilibrium of an element of volume. 

(a) General form of equations in polar coordinates. The polar co<3rdinates 
being r, 6, (f>, let the corresponding displacements be 

the cubical expansion A , and the three components of rotation nr^ , vx^ , rsr^ . 
Also let the components of total bodily force (per unit volume) be pF^, pF^, pF^. 
Then the equations of equilibrium for a volume element are the following : f 

5A 2u r 5 5t»r,"l 

(117) (X+2^)^--^-^[^(^,sin^)--^J=-,i^., 

1 aA 2u VdTiy d T 

(118) (X + 2M)-^g-^,[^--^^(...,sin^)J=-,i^„ 

1 aA 2u, r 5 dm '\ „ 

(119) ^^ + 2^).-sin"^a0-v[a.(''-«)-^J=-''^- 

* In the manner employed by Laemoe in tlie paper cited above. 

fSee Leva, Treatise on the Mathematical Theory of Elasticity, 2d edition (1906), p. 138, for the 
first members of these equations ; p. 56, for the values of A , VTr , We , •or^ . 
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while the cubical expansion and the components of rotation are given by the 
following equations : 

(121) 2^, = ^^^[^(™,sin^)-^(n*,)], 

(122) 2^, = ^-,[|^'-|,(«,sin^)], 

1 r 5 5m "1 

(123) 2.., = ^[^(..,)-^J- 

(6) Assumptions. The potential of the disturbing forces, being a spherical 
solid harmonic of degree i , is proportional to 

in which S^ is a spherical surface harmonic, and therefore satisfies the partial 
differential equation 

(124) sh!^^^ + si^|(-^^^) + ^(^ + l)'^. = 0- 

It will be assumed that 

(125) u^ = uS„ 

dS. 

(126) ug = v-^, 

1 dS, 



(127) '^* = ^Sn^a0' 

in which u and v are independent of 6 and 4>. These assumptions give by 
virtue of (124) 

(128) A.[li^-i(i + l)?]5„ 

(129) 2td-^ = 0, 

(ISO, -.4[»-^-^']a^.f. 
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and reduce (117), (118) and (119) to the following: 

(182) [,x+v)^,(^/-^>-.-(^+i):-) +i(.+i),'^(*i>-»)].,=-.^„ 

(133) [(^ + 2/.)^-,(-,J^^-<^+l)-j + -^^(^AJ-.jJ^-=-,i.,, 

(c) Bodily forces. We now proceed to express the values of (1) the dis- 
turbing forces, (2) the gravitational force due to the change of density of the 
attracted element, (3) the attraction due to the change of configuration of 
the attracting mass. 

(1) Disturbing forces. Let the potential of the disturbing forces be 

(135) Tf; = -|c,aJ^; 

then the terms in pF^, pFg, pF^ representing these forces are 

d W. r*-' 

r*-' d S. 





^ dr 






pdW, 
r 36 


= 


p 


8W. 





.-.OCX r —i 2 r*-' 1 3S, 

(^^^) r^^ ^d^ = - 3"^ ''P^" -^ sln^ # • 

(2) Attraction of original mass-configuration upon increment of mass of ele- 
ment. The attraction of the original mass-configuration per unit mass upon an 
element distant r from the center is 

r 

The increment of density is 

(139) -pA— p[^i<l^-i(.-+,):;]^,. 

Hence the attraction per unit volume, due to the increment of density of the 
attracted element, is 



(140) 



r V du u -.v^lr, 



1910] A GRAVITATING ELASTIC SPHERE 241 

(3) Attraction due to increment of density of attracting mass. The potential 
at a point (r, d, <f>) due to a shell of density D'S^, radius r and thickness dr is 

(141) ~^-j^^, D'S.dr (.>.'), 

(142) Jn^^^_.j)'S,dr' (r<r'). 

zi + 1 r 

On replacing D' S.hj 

the potential due to the change of density of the sphere is given by the expres- 



sion 



(143) 



(2i-f 1) a Lr'+ij„ \'('+^K' r'' dr' J"^ "^"^ 



r {•,■ ^s""' 1 d(r'^u')\ dr -[ 



To this must be added the potential due to the surface inequality of density p 



and thickness u, S. : 

1 t 



(144) 2nri-i^^- 

Adding (148) and (144), and reducing by integrating by parts the last term in 
each definite integral, we see that the total potential may be expressed in the 
form 

_ --(^»[,4,jr((^..)^.^)/-- 

The components of bodily force due to this potential are (per unit volume) 



d Y. 3w 



dr ( 2i + 
(146) 



r 36 
(147) 



+ ^(^+l)r-J"(^^;-^4)-J; + (2^ + l).]^, 
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(148) 

+ (z + i)r-] (v--7)^'-JsoV- 

(4) Components of total bodily force. Combining all the bodily forces, we 
find 

(150) 

while pi^^ differs from pF^ only in the substitution of l/8mff-BSJd4> for 
dSJdd. These values substituted in (132), (133) and (134) give the three 
equations of equilibrium for a volume element. It is seen that the factors 
involving and <l> cancel, also that the second and third equations are iden- 
tical. Hence the problem is reduced to the solution of two simultaneous ordi- 
nary differential equations for determining u and v in terms of r. 

(d) Final form of equations of interior equilibrium. We now get the equa- 
tions in final form by substituting in (132) and (133) the above values of pF^ 
and pFg . In writing the equations we shall replace u and v by new functions 

(151) e =: ~, a = -, 

and introduce as independent variable x = rfa . The equations then take the 
following forms : 

(152) +2i{i+l)fJM + pga\ — hx'+ (x-^ + 6e — i{i + l)a\x^ 



+^Tr(-^Xt(^+i)«'+^'>'''''^="'+^X^^«'-^'):c^-'^)]=^' 
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d& Q CL do, 

{\+fi)x2 — |-(3\+5//.)e+/xcc^ j^ +4iu.x~i- — li{i + l)\+2{v'+i—l)fi]a 

(153) +pcja^- |x'+ 2^- (^-^jT"^ [(^ + l)a'+e>"-^W 

§ 3. Surface conditions. 

The stress upon the surface of the actual strained body is assumed to be . 
This is equivalent to assuming the body to extend to or beyond the original 
boundary »• = a in every direction, and assuming the stress on this spherical 
surface to be a normal stress equal to the weight of the surface inequality of 
thickness u^S^ or e^aS■, e^ being the surface value of e. 

That this is correct to the first order of small quantities may be shown by the 
reasoning outlined above in the solution of the case i = 2 , which need not be 
here repeated. Assuming its correctness, the boundary conditions are the 
following : 

fr = - pgae^ S, 

when a; = 1 . 



du 



(154) 










rd = 
r(f) = 








In these 


we 


must 


put 


















fr 


= \A 


+ 


2/i 



(155) 



\ cr r r sin <j> d<pj 



r r sin ■ 
By means of (125), (126), (127) and (151) these take the forms, 

(da \dS, 

{ da \ 1 dS. 
\ ace J smff o(f) 



(156) 



so that the three equations (154) reduce to the two equations. 
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de 

(157) (X + 2/i)a!-T-+ (3X + 2/x)e — i(i + l)\a— —pgae, 

(158) ^J + « = <>' 

which must be satisfied for x = 1 . 

[It is easy to see what equations must replace these if different surface con- 
ditions are given ; for example, if there are known surface stresses proportional 
to S., dS./dO, 1 jsin -d S.ld<f), ov given surface displacements of a similar 
type.] 

§ 4. Solution of equations. 

A solution of (152) and (163) satisfying (167) and (158) is obtained by 
assuming 

(159) e=i:A.^. «=E^„a;'". 

The substitution of A^^x"^, B^x'^ for e and a yields in (152) the terms 
|[(\+2M)(mH3m)-Mi(i + l)]J^-J(J + l)[(X+/^)m-2/.]5^}a;"' 



(160) +W«([("^+6)^™-»(»+l)-5.]a'"+' 

+ 2^-l-l L~ m+J-1-3 + ^_i+2 * -m-i+2^j\' 

and in (153) the terms 

{ l(X+fji)m + (d\+5,i)-]A^+ lfi(m'+3m)-{i{i+l)X+2(;P+i-l),i)-]B^\x'^ 

(161) +^^''2i + lL m + i + S "" 

(i + l)(^ _j5 ) (i + l)(^ -j-B ) .-] 

m — z-f-2 m — i + 2 J 

Notice first that the terms in x' in these two expressions are in the same ratio 
as the original terms of degree i in (152) and (153) ; therefore, the terms of 
degree i in both equations can be made to vanish by a single relation among the 
coefficients in the assumed series (159). 

Again, for the terms of lowest degree, the coefficients of a;"* in (160) and (161) 
must vanish, and this is possible without the vanishing of A^ and B^ only if 
the determinant of the coefficients of A and B vanishes. The value of this 
determinant is found to be 

(162) (\ + 2fi) fi{m — i){m — i + 2){m + i + l){m + i + S) . 
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The four values of m for which this vanishes, with the corresponding values of 
B^l A , are as follows : 



m 


B^ {\+2,i){m^+Sm)-fii{i+l) 
A^ i(i + l)[(X + M)»i-2M] 


i 

i-2 

-i-S 


(i + 3)X+(i+5)ya 
(i + l)[a + (i-2)ya] 

1 

i 

(i-2)X+(i-4);a 

1 

i + 1 



Negative powers of x being excluded by physical considerations, the series may 
start with terms of degree i — 2, provided i is not less than 2 . If i is positive 
and less than 2, the lowest degree available would be i. 

We here assume i equal to or greater than 2, and let the series (159) take 
the forms 

e = ^;_2a;-' + A.x' + A+2^''^' + • ' '' 
(163) 1 

« = T A._^x'-' + B^x' + 5,+„«'+' + • • • . 

From (160) and (161) it is seen that (excluding the terms of degree i) the 
coefficients of a;'"+^ resulting from the substitution of (163) in (152) and (153) 
will vanish if the following equations are satisfied : 



[(X+2/.)(m + 2)(m + 5)-K^+l)]^.+2-KHl)[(X+M)(m+2)-2M]5„ 



(164) 



+ P9a 



'4.8m-(i-3)(i+4) 



[(m + 3)^^-KM-l)i?J = 0, 



(165) 



{m + i + S){m — i + 2) 
[(\ + ^)(m + 2) + 3\ + 5;.]^^^, 

+ |^(m + 2)(m+5)-[H^+l)X+2(^= + ^-l);^]|5, 
3 



m+2 



+ pga 



[m + i + ?,){m — i+ 2 



[(m + 3)^„-i(i + l)5J = 0. 



By means of these equations the coefficients A^^^, -^i+ti '■■■> -^,+2' ^.+4' "■ 
may be expressed in terms of A^, B^, leaving thus far arbitrary the three con- 
stants A._2, Af, B.. 

It remains only to equate to zero the coefficients of x% which as already 
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remarked requires but a single equation. The two differential equations will 
thus be satisfied by values of e and a of the assumed form, containing two 
arbitrary constants, provided the series thus determined prove to be convergent. 
The two arbitrary constants suffice for the satisfaction of the surface conditions. 
The solution of (164) and (165) gives 

A - _ C^^ 
"■+2 X + 2/it 

(166) 

[(m+4 )(m+6) M-^•(t+l)(X+2/.)][(m+3)^^-^(^•+l)i?J 

^ (w— i + 2)(m+i+5)(TO— i+4)(m + i+3)/tt ' 

(\R'l\ B Pga[m(X+^) + 5X + 4/x][(m + 3)^„-i(t + l)^,J 

l^-"*; ^m+2 (x+2ju)(m-t+2)(m+i + 6)(m-j + 4)(m + i + 3)A<.' 

(168) (-^+^l-^»+^^(^ + l)^».+2_ "^ + ^ W« 



(m + 3)^^-^(^ + l)5,„ („j, _ i + 2)(m + ^ + 3) X + 2/.- 

By means of these equations the values of all coefficients after the first two in 
each of the assumed series (163) may be expressed in terms of A^ and B^. 
These values may be expressed as follows, n taking values 0, 1, 2, • • • : 



(169) 



^.+2n+2=(-^r{Ki + l)^+2i(*+l)-Ci+2n + 4)(i+2« + 6)] 

(^ + 4)(^ + 6)...(j + 2« + 4)[^•(^+l)5.-(^ + 3)X] 



(170) 



-" (J+4)[2 • 4 . . . (2w+4)] [(2i+3X2i+5) • • •(2i+2ra+6)] ' 

^.._,o_..2=(-&)"+T(^+2«+5)-+^ + 2?^+4l 

V (i+ 4)(.-+6)- . ■(■i+2 >t+4)P( t+l)^,-(^+3)^ ,] __ 
^ (*+4)[2-4...(2/i+4)][(2i+3)(2i+5)--.(2i+2«+5)]" 



«+2n+2 



The coefficient of a;' may be found by substituting in (160) or (161) succes- 
sively w=i — 2,i,i + 2,---, remembering that iB^_^ = A^_^ , and that, as 
already pointed out, a single equation causes the vanishing of the coefficients of 
33* in both (152) and (153). This equation may be written in the following 
form, after dividing through by X + 2/i and writing h for gpaj{\ + 2/x) : 

*^'+ -^1^^+27^ "*■ i+1 ~ 2(2i+l)7^^ ~ V ^+2/^ 2(2^1))^'] 

+ 2i + l '-2- 27 + 1 ^ 2;rr4 ^&c = U. 

Also, from (169) and (170), 
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(172) 

[(^+4)( H6)---(i+ 2w+4)][(^+3)^-i(i+l).g,] ^ 

^ (i + 4) [2-4 . . • (2re+4)] [(2i+3)(2i+5) • • • (2l+2n+5)] ' 
by means of which (171) is reduced to a relation among the three constants 

The differential equations (152) and (153) are thus satisfied by (163), pro- 
vided the coefficients satisfy (169), (170), (171) and (172); and two of the 
coefficients still remain arbitrary. By properly assigning these, the surface con- 
ditions (157) and (158) may be satisfied. The solution may be expressed in 
fairly compact form by means of the following functions : 

^(^) = I:(-l)"[^(j + l)^+2i(^+l)-(i + 2w+4)(i+2n+6)l 
(i + 4)(i + 6)...(i + 2w + 4) 



X 



(t+4)[2•4..■(2;ft+4)][(2i + 3)(2i+5)•••(2^^-2?^-h5)] ' 

^(«)=E(-l)''^(» + 2« + 5)^ + ^ + 2r^ + 41 

(i + 4)(i + 6 )---( t + 2n + 4y ^ 

^(•i+4)[2-4...(2w+4)] [(2^ + 3)(2^+5).••(2i+2?^+5)]^ ' 

J^{z) = i:{-lY\i^ + Si + 2{n + S)\ 
L A' J 

(173) 

li±l)(i±^liL:(i± 2?L±±) n 

^(i+4)[2-4...(2n + 4)] [(2i+3)(2i+5) •• •(2^ + 2n•+5)] ^ ' 
L'{z)=J^{-iyli{i + l)'^+2i(i+l)-{i + 2n-\-4.){i + 2n+6)\ 

(^• + 2w + 2)[(^ + 4) (^ + 6).^■(^• + 2> ^ + 4)] ^ 

^(i+4)[2-4.-.(2w+4)][(2i+3)(2i+5)...(2i+2n+5)]^ ' 

(i + ^rij^-JKii + 4) ( ^ + 6)^^^(^ + 2/^ + 4)] ,, 

^(^+4)[2■4...(2?^+4)]"[(2] + 3)(2^+5).••(2i+2» + 5)]^ " 

The values of e and a, and the three equations for determining A._2, A^ and B^, 
may then be written as follows, the third equation being equivalent to (171), 
the fourth to (157), and the fifth to (158) : 
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(A) e = A._^x*'^ + A.x' + [(i + 3)^. - i{i + 1)^J hx'+^L{hx'), 

(B) a = v^i_2x'-2 + B,x' + [{i + S)A^-i{l + l)B^'\bx*+^M{bx^), 

xl-l + {i + l)bL{b) + {i + l)bM'{b)-\ = 0. 

These correspond to equations (A), (B), (C), (D), ( -E") obtained in Part 
I of this paper, and reduce to them when i = 2 . It must be noted, however, 
that the notation is slightly different, e, a and the functions (173) differing 
by constant factors from the corresponding functions in Part 1. 



